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ABSTRACT 


iminvwestenedtron 1s mide into the accuracy of the description of the 
actual mixing-ratio profile by a power-law approximation. The approxima- 
tion is obtained by a least-squares technique in order to derive the best- 
fit exponent A of the pressure profile. The <A -values were computed 
for soundings at 12 geographically diverse stations over the United States 
emerne, the period 16-23 March, 1971. The A -values at each een 
were found to undergo a time-variation based mainly on the synoptic-scale 
variations in the mixing-ratio profiles. A stepwise multiple regression 
procedure involving up to four variables from the temperature-humidity 
soundings was utilized in order to ''predict'"' the value of A _ from gross- 
parameters of the soundings by station. Tests were performed at four 
Stations in comparing the observed precipitable water vapor Wy with that 
computed from the power-profile Wes which depended upon A computed from 


the sounding. The values of W., gave a high correlation, R(W Wa emo? 


cs? 
for each of the four stations. However, the values of W.. using A from 
the predictive multiple regression equation gave somewhat larger standard 


errors at each of the four stations. The Wom7Values when compared to W, 


indicated some measure of predictive skill. 
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I. INTRODUCTION 


Reitan (1963) established a correlation factor of 0.98 between the 
logarithm of the monthly precipitable water vapor and monthly dewpoints 
based on 540 observations at numerous locations. The regression formula 
used by Reitan was 

inigieeee + Bt e (1) 
where A = -0.981 and B = 0.0341 are the best fit coefficients resulting 
from the regression procedure. 

Smith (1966) has shown that an equation of form (1) may be derived 
from the well known empirical equation of Tetens (1930) 

e =E, x 19 (4 ty - B)/(tg + ¥ ) ©) 
provided that the mixing-ratio profile may be parameterized. The constants 
are Eo = 6108 dynes cm7™“, a = 7.5, B = 238.1 F, and yY = 395.1 F. The 
required parameterization of the mixing-ratio can be accomplished with 
reasonable accuracy by the power law approximation 

w = Wo (p/Po)” (3) 
By combining eaaciione (2) and (3) with 

Wo = 8 ( &5/Po) (4) 
Smith was able Be show that 

Pee omuie ete cane 2 026) (at, = By/440 =... 

(5) 

-...-ln( A+ 1) 

Since climatological values of the exponent A are variable with latitude, 
Smith concluded that the coefficient A of (1) is likewise dependent upon 


Mati tude in Reitan's formulation. 
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TABLE I 


A -values from Smith (1966) 





LATITUDINAL ANNUAL 
ZONE -ON WINTER SPRING SUMMER FALL AVERAGE 
@ - 10 Baan 2.85 2.80 2.64 2.91 
t= 20 2.99 BeOg 2.70 2.93 al 
Zar= 30 3.60 3.00 2.98 2.93 Bail 
30 - 40 3.04 Baill 2.92 2.94 3.00 
40 - 50 2.70 2.95 2h Deal 2.78 
50 - 60 E52: 3.07 2.67 2.93 2.79 
60 - 70 1.76 2.69 2.61 2.61 ae 
70 = 80 1.60 1.67 2.24 2.63 2.03 
80 - 90 eat 1.44 1.94 Defy) 1.62 

NORTHERN 

HEMISPHERE Day 2.64 2.62 DT 2.61 

AVERAGE 


iy 





Bolsenga has attempted to adapt Reitan's formulation (1) to daily 
and even to hourly observations and finds that the explained variance of 
total water vapor falls off to 72% and 64% in the respective cases. This 
could have been anticipated since total water vapor depends in each case 
upon the entire mixing-ratio profile, which depends less upon surface 
dewpoint when a decreased time interval of sampling is employed. 

The equivalence of equations (1) and (2) to the A -profile (3) holds 
best, according to Smith, when applied to climatological data. Climato- 
logical values of XA were calculated by Smith (1966), from London's global 
mean radiosonde data (1957), with resultant A -values presented as a 
function of latitude and season as shown in Table lI. 

The concept of an accurate mixing-ratio moisture profile from sound- 
ings composed of few or widely spaced levels is a very valuable tool. 

It is especially important in the field of radiation flux measurement 
where values for the mixing-ratio (moisture profile description) are 
needed at all levels to the top of the atmosphere but are frequently not 
available from sounding levels much above 500 mb. Recently Plante (1973) 
and Martin (1973) have extended the use of Smith's power-law profile to 
one of the form 

Ww = w599(p/500) * (6) 
with the reference or key level at 500 mb rather than a surface reference 
level which has a diurnally variable mixing ratio. Martin (1973) using 
the Fleet Numerical Weather Central, Monterey, California, computer print- 
out of mixing ratio at six pressure levels in the vertical at FNWC erid 
points found a correlation coefficient of .998 or higher between observed 
mixing-ratio and a mixing-ratio computed by Eq. (6), based on A -values 


calculated at each grid point for the 75 grid points considered. Martin 


ne. 





concluded that the power-fit as described by Eq. (6) is accurate especially 
with soundings having widely separated reporting levels in the vertical. 

In this thesis the following questions will be investigated: for a 
given station for a given period of time, can a power-fit of the nature 
of (6) be obtained which is representative of the actual moisture profile 
at 50 mb-level increments, and what degree of accuracy can be expected? 
Further, the conditions which would tend to invalidate the accuracy of 
the method will be delineated. Two methods of verification of the results 
will be considered: 

(1) the explained variance of mixing-ratio resulting from the use of 
power-law profile, and 

(2) the accuracy for the calculated precipitable water obtained from 


the profile-parameterization compared to that observed. 


ES: 





It, DATA PROCESSING 


A. THE ORIGINAL DATA 


For the water vapor analysis, twelve stations were selected on the 
basis of their locations and the availability of data. The time chosen 
was an eight day period from 16 March to 24 March, 1971. The station 
locations, elevations, and mean surface pressures are shown in Table II. 
As can be seen from Figure 1, the stations selected lie in a band across 
the southern, and in a second band across the northern continental United 
States. The data was obtained from NWSED, Ashville, North Carolina. All 
data levels were selected at levels corresponding to integral multiples 
of 50 mb (up to 400 mb), in addition to the surface level. 

The following analysis was used as the basis of solving for A , 
Biareime from Eq. (6): 

tos utp) = HOS see (7) 
Here p ranged from the surface value to p = 400 mb, and corresponding 
values w(p) were computed at each level p using sounding-level data. If 
the ieft side of Eq. (7) is denoted as Y, then this equation becomes 

Y= AX | (8) 
where X = log (p/500). In the average sounding which extends to 
p = 1000 mb, there would be 14 simultaneous values (Y;, Xj) so that only 
a best-fit solution for A is feasible. 

This amounts to solving the linear regression equation, Eq. (3)e0 Lou 
the coefficient A of the independent variable, with the condition of the 


zero-intercept forced on the constant term. In the program library of the 
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W. R. Church Computer Center, the most convenient program for performing 
this objective was BIMED 02R (Dixon, 1966). While this program is primarily 
Mised for Stepwise regression in the case of multivariate regression analysis, 
it also has the desirable option of making the required transformations of 
variables first to ratios, and then to their logs as in Eq. (7). This 
meguires that the sounding data has been prepared in properly formatted 


form on IBM cards. 


B. CALCULATION OF MIXING-RATIOS AND OF A 


Mevine decided on the method of handling the regression, the next 
computation was that of mixing-ratio to be specified by the power-profile 
mee, Eq. (6). 

The data as it was received (by courtesy of Mr. Russell Schwanz), was 
card-punched in a form containing the pressure, temperature (C), and 
relative humidity at 50 mb-intervals from the surface to the top of the. 
moisture sounding (which was generally at 400 mb). The next step then 
was to compute the mixing-ratio in its approximate form 

w = .62197 (RH) e. (AL), (9) 
meme €, is the saturation vapor ere computed over a plane water 
surface for T 2 OC, and e, 1S computed with respect to equilibrium 
over a pure, plane ice surface at temperatures T< 0C. In doing this at 
each sounding level, the two versions of the Goff-Gratch equations for e€sy 
and esi, respectively, were used in the forms indicated below as Eq. (10) 
mad Eq. (11). 

logi0 esw = -7.90298 (Ts -1) + 5.02808 logig (Ts)- ... (10) 

es co 

Popes com On?) (1OLL.334(1-T/Ts) -1) +... 


ces 2 oes (l0.2- 49149 (1./T-1)-1) + lopi9 ews 
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Logyoes; = ~-9-09718(To/T-1)-3.56654 logy9(To/T) + ... (11) 
meOrer cor oo (red, I) + logy eyo 

where T, = 373.16 K and Lt Los o KK. 
In (10) esw is the saturation vapor-pressure of pure nS the steam- 
point temperature and esi is that of pure water at the frost-point 
temperature (IT <0 C). Then w(p) is evaluated from either equation (10) 
or (11) utilizing Eq. (9). Once the values for w were obtained at each 
level in each sounding the best-fit value for A was obtained in the 
manner suggested in connection with Eq. (8) and Section III (A). Thus a 
unique value of A resulted for each sounding. The A -values for each 
Station were then grouped with other synoptic sounding data for each 
individual sounding in order to account for the reduced variance of A 
for the sounding-sample for the particular station. In this connection, 
the stepwise regression program BIMED 02R was used to determine the 
explained variance of A at each station. Comparisons were made of 
computed w-values resulting from Eq. (6) versus observed w-values over 


the sounding station. 


C. CALCULATION OF PRECIPITABLE WATER VAPOR 


Precipitable water was obtained by integration for both the observed 
soundings and ore the corresponding <A -profile cases and a comparison 
was made for several stations. The formulas used in the calculations are 
shown as follows. 


a. The observed value of precipitable water vapor is 


Pia 
W = “dp = 1 put A Pb ig cy ee (12) 
g 980 2 eS 2 DES 
Fi 
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Z 25 2 1s es. 
b. The model value of precipitable water vapor is 
Ps Be 
r 1 _ 
a “SOE, ba) £ dp = W500 (500) ( p yn 
500 GA ol) 500 
Prop : : 400 


(13) 


W500 (500) Po\ A +1- “0 yen 
980 (A +1){ \500 500 


The comparison of the sample sets of the two precipitable water values 
will be used as a second measure to determine the accuracy of the A - 


profile method. 


ie, 





Pi eeolatisStiCAL CONCEPTS 


A, THE LEAST SQUARES FIT 


In this thesis, the least-squares fit was employed in at least two 
contexts. First, it was employed in solving Eq. (8) for XA at each 
Station and sounding time. Here the variables Y and X are given point- 
wise by 


X, = log _F and Mera Loe we OS ills) 
500 





W500 
with i representing a level numbered from one at the surface, to two at 
the next higher level (where P first becomes an integral multiple of 
50 mb), etc. Standard texts in statistics (e. g. Brownlee, 1960) show 
that the parameter A of the power-profile as in Eq. (8) can be deduced 


from the best-fit condition 


>> X;? (16) 
i 
If the zero-intercept. had not been required in (16), A would have been 


estimated from 


eae Oa = Y) 


eh 
N Ox (17) 





and the constant term of the best-fit simple regression would then be 


Dey: 


N (18) 
where N is the number of sounding-level observations entering into the 


power-profile (normally here N = 14). 
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The combination of conditions (17) and (18) gives the standard best- 
fit solution for the simple independent variable case 
Veen A OX cio) 
Equation (16) gives the corresponding result when the zero-intercept 
(Gq X; = 0) is required. In either case A is related to the simple 
correlation coefficient R)72 and the variables Y and X through the standard 


formula 


Ri2 = ox (20) 


Hence a value of Ryjy close to unity means an effective least-squares fit 
by the power-profile and simultaneously a reliable value of A. In 


Eq. (20), the statistical parameters 9. and © y are defined as 


aie 2 
2 o. = ee os)! 
N 
C2) 
b 0, = 


au _ sf) 
N 

Nete that %, and Oy of (21) are applicable to the zero -intercept case 

and apply to the more general case of the standard regression formula 
re Si R12 _7Y (x; - x) (22) 

wk 
eigen! Y; is replaced by Y; - Y and Xx, isi replaced by X; - Ke 
It was possible to apply BIMED O2R regression procedure to the 

14-values of Xj and Y; per sounding (that is, with each sounding listed 

on 14 cards) and to generate not only the A -value for that sounding, 

but also the simple correlation R12, which is a measure of the "goodness 


of fit" of the power-profile for the particular sounding case. Such 
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results for A and Rj? have been tabulated under each station and for 


each observation time (Tables III through XIV). 


B. THE STEPWISE REGRESSION PROCEDURE 


If there are N samples of the statistically related variables 
4a 34> ---aX,a) i =1,...,N, it is possible to apply the least- 


Squares concept to minimize the least-squares form 
3 Z 
oT A = Ay = AgX9q Ee lee Seni) = minimum C23) 


i=l 
The procedure of selecting the least-squares fit of coefficients to 
minimize Eq. (23) is called multiple linear regression. In this process, 
there is no assurance that the variables Xo, eee Xy used to explain the 
variation of Y, are the most effective, and, in fact, some of the X; may 
meeecontribute to the explained variance in any significant manner. 

BIMED O2R is a statistical procedure for developing a multiple 
regression equation in which the variable x. added at each step accomplishes 
the greatest reduction in the unexplained sum-of-squares, Rd)". If 
k-1 variables are to be added in a stepwise manner the first one added 
is that which has che lessee Sai plem corre latcGron COCLLICLent Rij ot all 
k-1 possible comparisons. As a result, the first estimate to the desired 
multiple eepytaRS LOR equation may be obtained in the form: 

Cul teen I) (24) 
feeeimag7, 2), selected as in (17), (18). A property of this first-step 
selection is that the sum-of-squares explained by Eq. (24) is given by 


; (S.S. Expl) ] = RGD | on-Fi5+..-+0ran 4)" | (25) 


A result which is proved in standard statistics texts (e.g. Crow et al, 


960) . 
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The BIMED OZR program then forms the data-sample residuals. After one 


step, all values have the form 


7 1 
"OA a Ge ee Gt (1) (26) 


The program then examines all of the simple correlations of uP ea rae 
X, (the remaining independent variables) and selects that which has the 
Paenest Simple correlation coefficient among RW 4X), ies Oana en 
Since this simple correlation is based upon the residual after the effect 
fumveriable X5 has been removed (see Eq. (26) ), it is typically called 
Mreminrst Order partial correlation with the notation Rj3.2, the under- 
standing being that the variable having this highest correlation at 
step 2 and thus being selected is called X3. Then the regression equation 
relating ee X3 is obtained as described in connection with Eqs. (21), 
(22), (23) and assumes the form 

mee by, + by3 X3 (27) 
As a result of the eee simple regression, Eq. (27), the sum-of- 


1 
squares of Yy (not yt ), 1s now Statistically explained in the amount 


Pi Bee =e 


The statistical parameter Rj 93 is the so called Mmiwenp be cornelation at 
step 2, that is to say of variableY, depending jointly on the variations 
of both X» and X3. These independent variables were selected according 
to the criteria of having the largest simple correlation, and the largest 
/meest-Order partial correlations respectively. 

itesmultiple correlation coeitficient satisfies the important chain- 


rule of multiplication 


2 2 2 
E> By 23 = (TRy>) CR,3 9) (29) 
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and after k-l steps, this rule assumes the form 


i= eee = (1 - Ryo) C= See - ‘Sie ee (30) 
Just as eee describes the fractional explained sum-of-squares of Y in 
Eq. (28), so ee describes the fractional unexplained sum-of-squares 
about the regression surface. This parameter has decreased relative to 
mice Ok step one, provided ee > Ryo. The standard error of estimate 


is related to the multiple correlation coefficient by 


2 2 2 
Syix = Wy (1-Rj2) Cr) 


The parameter standard error, Sy x: Wimelebpe listed in ftapres Tll..4.,XK1V. 
aimee), the subscript Y!X denotes the best fit of Y on all X,. 

The most attractive feature of the stepwise regression is embodied in 
the chain-rule of Eqs. (29), (30), with each step making the selection 
Of the highest partial correlation coefficient for use in Eq. (30). The 
result is that when the unexplained fractional sum-of-squares does not 
become significantly smaller, the process of adding variables to the right- 
hand side of 

2] = Ay a0 Ayko a A3 X3 agree A) Xy, (32) 
Should be terminated i the (k-1)th step. 

In Section IV, application of BIMED 02R will be used on the 16-sample 
values of A derived at each individual station. This is done since it 
turns out that a time variation in A occurs with the synoptic situation 
at individual stations. In some cases this variation, as indicated by 
the standard deviation @, , is so small as to be limited by the accuracy 
of the data. In other cases, a rather large variation of A occurs and 
meewas desired to comsider how the time-varying structure of the soundings 


contribute to the A -variation. For this further study each X -value 
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was added on an IBM card containing certain relevent parameters of the 
sounding from which A was derived. For all soundings and stations a 
set of six independent parameters (for example T(850) and w(850) and 
others to be described in Section IV) were included along with the A - 
value for the identical sounding. 

Thus BIMED 02R was also usedinthe important stepwise sense in 
describing the <A -variation in the sounding structure over the set of 
the 16 soundings per station. The detailed discussion of the variables 
used in the regression for A is given in Section IV (B, C), together 


meee the discussion of the results. 


Ge ESTIMATES OF STATISTICAL SIGNIFICANCE 


l. Of A _-Profiles at Sounding Stations 


Here the correlation coefficient Rio of the regression equation 
(7) is the basis of the estimation of the accuracy of the A -profile. 
The significance of the test is based upon Eq. (25) which gives the 
result 

: ; . 2 
Fractional Explained Variance = Rj5 (33) 
whereas the fraction residual (or unexplained) variance is given by 
2 


Fractional Unexplained Variance = 1-R), (34) 


The significance test is based upon the F-statistic (which is listed in 


the BIMED O02R output) 9 
R 
rE 13) = 2 (35) 
(1-RZ yikes 


The critical value of F(1,13) at the 99% confidence level is listed in 


Pedttastical tables as 
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He ee, =O] 


This 99% confidence level is satisfied whenever Rj) 5 of Eq. (34) exceeds 
Rio (eritieal ~) 2 641 


ine) but one of the A -profiles, Rio exceeded this value, and in this 


exception R} 9 was 0.614, corresponding to a confidence level of 97.5%. 


2 
In summary, the statistic Rj) listed for each sounding at the 12 


Stations, is the explained fraction of the variance of Y in that sounding 


iseemrables III,..., XIV). 


Finally, BIMED O2R lists also the standard error of the coefficient 


Z 
Seandard Error Of A. = 9A / TR - N= 14 (36) 
Z 
Rj > (N-1) 


This statistic appears in column 3 of Tables III,... XIV, and describes 


the sharpness of the heer te i. C7) 4 


2. Sienificance Tests of the Temporal Variation of A 





Here the fractional explained variance may be computed from 
Eq. (28) after the kth entry of an independent variable - drawn from the 


BEyuopeic SOUNding variables has been included in 
ees A Ka eet Al Ky (37) 
The fractional explained variance is given by R* in accordance with 


(28), whereas the fractional residual variance is 1-R*. In general when 


a9(A) is sizeable, it is possible to generate a multiple regression equation 
of form (37), which significantly reduces the residual variance. In 


general the criterion used for entering additional variables for explaining 


the  X -variance was based on a stepwise selection which added at least 


2% to the total explained fractional variance. This strategy was used 


meeause the F=statistic for this case is 
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Re > Reed 
lee GL N15 10) as eee ema (38) 
Cee) (N= k = 1) 
and here the critical value of ee is difficult to surpass with N = 16 
observations. The added fractional explained variance at step k is 


identical to the numerator of Eq. (38), and is required to be at least 


Beeeat each step. 


rat 





eek oULLS 


A. DISCUSSION OF A_ -PROFILES 


This particular phase of the research began with computer cards 
containing the following individual sounding data: pressure (mb), tempera- 
ture (K), and mixing ratio (gm/kg). There were 16 soundings of up to 14 
levels at each of the 12 stations (Table ITI). 

The BIMED O2R program was utilized in a least-Squares context to 
solve Eqs. 7 and 16 for a best fit A for each sounding at each station. 
The "zero-intercept" was chosen in the correlation so that the computer 
Solutton to Eq. (7) was of the form of Eq. (8) instead of Eq. (19). A 
discussion of the statistical methods and parameters is given in 
Section III (A). Table III lists the best-fit <A -value for each avail- 
able sounding at Waycross, Ga., over the test period. 

Zt should be noted that the least-square procedure for determining 


from 


ee CG?) 


minimizes the sum-of-squares of the residuals > _ Hences the estimator 
Y may be defined as 


Y = A log P/500 = log w,(P) (40) 





W500 
where the observed ws5gQ appears in the denominator of Eq. (40) by the zero- 
intercept concept. On the other hand, the observed Y is defined as 


Y = log w/ W590 so that the residual becomes 


A 
¢, = Y - Y = log w,(P) (41) 
w(P) 
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In Table III, successive values of A obtained by A -profile fit 
[Eq. (39)] have been listed along with the standard error of each A 
determination. The small standard error of A indicates the close fit 
of the A values by Eq. (16). 

fmmeacdition, the fractional explained variance RY for each sounding 
is listed alongside R,9 under column 4. Finally, the standard error of 
estimate is listed in the final tabular column as defined by Eq. (31). 
The value of Sy;y is related to the standard deviation oy through 

Saicecl > Rio) oy (42) 
These ratios of Eq. (42) obtained by applying the values of Syjx and 
meomecolumns 5 and 6 of Tables III through XIV are generally of the 
Order of 0.3 or lower. There are exceptions, but these can be identified 
by the values of Rj» (generally Rio 2 .95) for the respective soundings. 


ime iable Iii, the standard error may be written in terms of the mean 


residuals - Squared as 
N , a 
Sy ix = 2 meone'c | 7 eaeealirwaeatian (43) 
aa W 1: E3 W 
where the final term in Eq. (43) represents the mean scatter about the 
A 


graph of Y= A X. 
For example, at Q0OQO GMT, 21 March, 1971, the standard error at Way- 


cross, Ga., was 0.1131. Therefore, .1131 = log w+ |AWol , hence it 


W 


follows that w+ |[AWcl = 1.298. This means that the mean residual ratio, 
W 
after regression, for this case is [|A%cl = 30%. Analysis of the actual 


residuals for the soundings shown indicates that the majormey Gf this €rren- 
ratio occurs at the highest altitudes recorded in the sounding where w is 


the smallest and where actual instrument accuracy is doubtful. 
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TABLE ILL 


Values of A and related statistics 
for the data-sample for Waycross, Ga, 





STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE A -VALUE 2 
DAY-GMT CALCULATED A-VALUE R12 R12 oy SY|X 
16-0000 2.9414 Or 555 953 .908 0.5861 0.1837 
16-1200* 1.5298 0.4277 .704 .496 .4130 .3043 
17-0000 4.5798 0.1950 .988  .976 . 8809 ee 
17-1200 2.7855 0.1246 987 .974 De TS .0888 
18-0000 22 OF TL87 .983  .966 577 0846 
18-1200 2. 6057 Orage .983 .966 BS05i7 .0978 
19-0000 2 3259 0.3495 none 2703 .5038 .2491 
19-1200 1.8176 0.2879 868 .755 .3981 . 2049 
20-0000 3.6279 0.2036 .980 .960 . 7030 .1447 
20-1200 2.5391 0.2308 950 .903 .5083 . 1643 
21-0000 3.08387 0.1589 .983 .966 2507 6 ELISA 
21-1200 Pea 250 0.1087 .987 .974 .4680 NOS 
22-0000* 0.9815 On3252 oe asl .2911 . 2316 
22-1200 3.2945 0.2850 $955, 4.902 ean 2 gal 
23-0000 4.3912 0.2107 985  .970 8473 .1499 
23-1200 4.5008 02295 982 .964 sor .1724 
TIME MEAN 228557 . 869 
Sap. DEV. 1.0401 
CA) 





(@erresponds to Ryo < .85. 








30 











TABLE Ly 


Valmeswof | A and related statistics 
Bor tne data-sample for Charleston, S. C. 





STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE A -VALUE 9 
DAY-GMT CALCULATED A -VALUE ES Sol oY SY |X 
16-0000 2.3199 0.1174 .984 .968 0.4485 0.0863 
16-1200 CS 7 valk 0.2933 .972 945 0.8494 0.2088 
17-0000 3.9920 0.1322 .993 .986 0.7646 0.0941 
17-1200 3.8305 0.1629 .989 .978 0.7379 0.1160 
18-0000 2.4849 DL De 970 .941 0.4883 0.1229 
18-1200 2.7495 0.1236 987 .974 0.5319 0.0883 
19-0000 2.9405 0.1308 .988  .976 0.5684 0.0934 
19-1200* 1.5187 0.4050 21 520 0.4013 0.2886 
20-0000 4.5054 0.2364 .983 .966 0.8711 0.1680 
20-1200 7575 0.2153 961 .924 0.5369 0.1532 
21-0090 4.7245 0.1834 .990 .980 0.9077 0.1305 
21-1200 2.4194 0.1502 976 .953 0.4727 0.1071 
22-0000 MISSING DATA 
22-1200 2.5801 0.2807 .931 .867 0.5281 0.2001 
23-0000 4.4411 0.2290 .983 .966 0.8597 0.1631 
23-1200 3.5360 0.1560 .988 .976 0.6809 0.1110 
TIME MEAN 3.2729 928 
STD. DEV. 0.9832 
mein.) 


*Corresponds to Ry» < .85. 
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TABLE V 


Values of A and related statistics 
for the data-sample for Lake Charles, La. 








STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE A -VALUE 2 
DAY-GMT CALCULATED A-VALUE R12 R12 oy SY 1X 
16-0000 D ONE 0.1461 .984 .968 0.5653 0.1041 
16-1200 4.0954 0.1065 996 .992 e7e37 .0759 
17-0000 4.5898 0.1810 990 .980 8835 1291 
17-1200 3.5190 0.2167 976 .953 6878 1547 
18-0000 4.0538 0.1916 986 .972 7843 1367 
18-1200 3.8884 0.1964 984  .968 7528 1399 
19-0000 2.7501 0.0960 .992 .984 5273 0684 
19-1200 3.0506 10.2307 965 .931 6030 1646 
20-0000 2.4835 0.0959 .990 .980 4782 0684 
20-1200 28351 0.1944 971 .943 5575 1388 
21-0000% esas 0.2333 .813 .661 2757 1665 
Bi 1200 1.6059 0.2349 "8¢5mee. 783 3462 1676 
22-0000 4.3386 0.1766 .992 .984 Tan 1041 
22-1200 5.2771 0.1789 .994 .988 8705 1016 
23-0000 4.4083 0.2264 983 .966 8532 1612 
23-1200 4.5449 0.3695 963  .927 8417 ens 
TIME MEAN: 34711 936 
STD. DEV. 1.1352 
g( A ) 





Beorresponds to Ri5 < .85. 
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TABLE VI 


Values of A and related statistics 
for the data-sample for El Paso, Texas. 











STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE A -VALUE 2 
DAY-CMT CALCULATED A -vALuE *12 %12 Oy Syitx 
16-0000 a Sy 0.3078 .964 .929 0.5391 0.1509 
16-1200 3.3290 meessa  o76 4953 0.5050 ue, 
17-0000 enc%6 0.1695 .988 .976 0.5195 0.0832 
17-1200 3.2455 aees47 986 2876 0.5123 0.1886 
18-0000 3.7536 0.2578 .997 .994 0.5658 0.1260 
18-1200 2 390 0.1721 .987  .970 0.4938 0.0843 
19-0000 3.3597 Gein cumeeoesn moo0 0.4990 0.0513 
19-1200 2.4246 0.2875 .936 876 0.3841 0.1414 
20-0000 2.5997 0.2698 .950 .902 0.4057 0.1327 
20-1200* 0.6696 0.1802 .778 .605 0.1311 0.0868 
21-0000 2.6959 ApeinlG nos9 ‘Vae2 2032 Ge 1524 
21-1200 3.1878 0.1860 .983 .966 0.4787 0.0911 
22-0000 We 213 0.2029 .982 .964 0.4995 0.0994 
22-1200 2.4686 aeosii) aoe 7342 0.4003 0.1722 
23-0000 2.3344 oeigs5 morse 4947 0.3545 0.0860 
23-1200% 0.7583 0.4089 .526 .277 0.2192 0.1965 
TIME MEAN 2.7771 871 


SDs DEV. 0.9147 
o7-(C A) 
“Corresponds to Ryo < .85. 
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TABLE VII 


Values of XA and related statistics 
for the data-sample for San Diego, Calif. 








STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE AX -VALUE 2 
DAY-GMI CALCULATED \ -VALUE Rig Ry ay SYIX 
16-0000 Boley .0897 .996 .992 segs .0637 
16-1200 3.8860 BLAa? 991 .982 7444 . 1024 
17-0000 3.4676 .0874 996 .992 . 6606 .0620 
17-1200 3.2899 AGT .988  .976 (5932 .0943 
18-0000 3.9886 Sone 987 .974 W205 . 1230 
18-1200 Olea a DE 992 .984 . 7684 los 
19-0000 3.6748 .1190 993 .986 7023 .0845 
19-1200 2.4739 70 OS qeaeisl 4817 Paislells 
20-0000* 1.9672 .3994 7o07 9 2650 4627 . 2836 
20-1200 3.9124 .1845 .986 .972 .7529 . 1309 
21-0000 33273 . 2063 976 .953 .6470 .1465 
21-1200 3.8186 . 1946 .984 .968 a7O mesa? 
22-0000 2.6873 2288 oe «| CON 80 . 1626 
22-1200 2.3642 23 .938  .880 .4786 nial 
23-0000 2.1203 sp Salk 919.845 .4378 .1790 
23-1200 2.8397 . 1020 992 984 .5434 n0725 
TIME MEAN 3.2096 .938 


oD, DEV. 0. 70i2 


g,( A) 


Begrresponds to Rig < .835. 





TABLE VIII 


Values of A and related statistics 
for the data-sample for Caribou, Maine. 











STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE A -VALUE 2 

DAY-GMT CALCULATED A-VALUE Ri2 Ry Oy Sy 1X 
16-0000* 1.3979 0.2714 .841  .707 0.2937 0.1661 
16-1200 4.3499 0.0839 998  .996 0.7705 0.0534 
17-0000 3.5403 0.1644 .988  .976 0.6588 0.1047 
ie7= 200 2.9646 0.1566 .984 .968 0.5350 0.1003 
18-0000 3.1736 0.2673 960 .922 0.5869 OMe: 
18-1200* e775 0.4351 P7238 523 0.3878 OnE 
19-0000 4.6004 0.1396 .995  .990 0.8453 0.0884 
19-1200 3.7880 0.2120 .982 .964 0.6859 0.1359 
20-0000* 1.0878 0.3522 .666 444 0.2901 0.2254 
20-1200 2.7297 0.2756 .957 .916 0.4660 Gra 
21-0000 5.0908 0.3520 .979  .958 0.9154 0.1960 
21-1200 51735 0) Bae .989 .978 0.9479 0.1461 
22-0000 3.7348 0.0533 999 .998 0.6789 0.0335 
22-1200 4.5463 0.1532 994 .988 0.8320 0.0965 
23-0000 6.5325 0.5259 975 9x95 1 1.1656 0.2745 
23-1200 4.5879 0.1398 996 .992 0.8014 0.0769 
TIME MEAN 3.6797 892 
STD. DEV. 1.4896 
Pe( A) 





*Corresponds to Ryo < .85. 
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TABLE IX 


Values of A and related statistics 
for the data-sample for Maniwaki, Que. 








STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE A -VALUE 2 

DAY-GMT CALCULATED A-VALUE Ri2 = Ry2 Oy Sy 1X 
16-0000 5.5686 0.3228 .981 .962 1.0019 0.2053 
16-1200 3.0656 0.2679 971 .943 0.4894 0.1246 
17-0000 3.6455 0.1713 987 .974 0.6561 0.1097 
71200 3.2185 0.2048 977.955 0.5867 0.1315 
18-0000 3.2674 @22283 976 .953 0.6138 0.1405 
18-1200 2.8017 0.1775 977.955 0.5597 0.1249 
19-0000 4.2382 0.1084 .996 .992 0.8075 0.0770 
ie 200% 1.1753 0.3235 .724 = 524 0.2892 On2077 
20-0000 1.8514 0.4214 .874 .764 0.3945 0.2075 
20-1200 3.0911 0.1692 .987 974 0.5533 0.0945 
21-0000 ones 0.0641 999 .998 0.5703 0.0320 
21-1200 eee? 0.1215 991 982 0.5029 0.0714 
22-0000 2.3573 0.2628 a33° 670 0.4481 0.1680 
22-1200 2. 6480 0.1180 .989  .978 0.4888 0.0747 
23-0000 3.4791 0.1644 .989  .978 0.6695 0.1038 
23-1200 3.8057 0.3301 971  .943 0.6607 0.1670 
TIME MEAN 3.1407 922 
STD. DEV. 0.9849 
7rC A) 


*Corresponds to Rio < ToD 
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DATE 


DAY -GMT 


16-0000 


16-1200 


17-0000 


17-1200 


18-0000 


18-1200 


19-0000 


m1 200 


20-0000 


20-1200 


21-0000 


21-1200 


22-0000 


22-1200 


23-0000 


23-1200 


TIME MEAN 


STD. DEV. 


aC A) 





Values of Xr 


X -VALUE 


CALCULATED 


2.4004 
2.4220 
4.3460 
3.6100 
geri 
3.5366 
2.0358 
2.2482 
2.8100 
3.2512 
33635 
Dales 
ORY, 
3.3318 
3.8431 


3.7412 


5. L6G 


0.6940 


STANDARD 
ERROR OF 


A -VALUE 


0 


750 


Pes: 


Syllve ils 


woo? 


ples 3 


WS Te) 


2896 


LILES, 


Jaye: 


0385 


1394 


2288 


~ 1823 


sgl ere: 


PEOLS 


Oy 7 


TABLE X 


and related statistics 
for the data-sample for Sault-Ste-Marie, Mich. 


CORRELATION 
COEFFICIENT 

2 
Rig = Ry 
.932 .869 
986 .972 
995  .990 
995 .990 
994 .988 
994 .988 
904 .817 
.943 889 
979  .958 
.999 998 
991 982 
nas. 893 
981 .962 
993 .986 
997 994 
996 .992 

954 


Bu 


STANDARD 
DEV IAT ION 


Cy 


0.4486 
0.4355 
Ore77 310 
0.6639 
0.6828 
O14 Gals, 
0.3889 
0.4225 
Ore ZS 
Ono71/ 
Ono221 
Oma / 
0.5614 
0.6384 
0.7333 


0.7496 


STANDARD 
ERROR 


SY1X 


0.1741 
Onur 70 
0.0780 
0.0697 
0.0807 
0.0721 
Oh 33 
0.1463 
0.1108 
0.0242 
Oe0S77 
0.1446 
Os Las: 
O27 72 
0.0640 


0.0680 





TABLE XI 


Values of A and related statistics 
for the data-sample for Bismarck, N. D. 








STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE A -VALUE 2 
DAY-GMT CALCULATED A-vaALUE 12 812 Oy Sy 1x 
16-0000 2.6549 0.2499 foo5e) | 59h? 0.5027 0.1564 
16-1200 3.4759 0.2024 984  .968 0.6213 0.1180 
17-0000 6.5097 0.1724 996 .992 1.1819 0.1080 
17-1200 NO DATA AVAILABLE 
18-0000 NO DATA AVAILABLE 
18-1200 3.6856 0.2127 984 .968 0.6601 0.1243 
19-0000 4.4516 Gn 78 986 .972 0.8179 0.1427 
19-1200 4.5117 0.1570 995 990 0.8724 0.0955 
20-0000 3.4355 0.1859 .983 .966 0.6154 0.1180 
20-1200 2.0076 0.2615 .918 .843 0.3888 0.1611 
21-0000 - 4.9089 Gmilozs 992 .984 0.8086 0.1088 
21-1200 4.5735 ° 0.1341 996 .992 0.7678 0.0743 
22-0000 5.0893 0.2730 987 .974 0.9794 0.1640 
22-1200 3.8804 0.0609 999 998 0.7341 0.0382 
23-0000 3.6326 Gae285 .954  .910 0.6710 0.2088 
23-1200 2.1195 0.1005 .987  .974 0.3789 0.0639 
TIME MEAN 3.9240 .960 
STD. DEV. 1.2147 
o.C A) 
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Values of A and related statistics 
for the data-sample for Glasgow, Mont. 











STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE X -VALUE 2 
DAY-GMT CALCULATED A -VALUE Rig =: Ry Oy Sy1X 
16-0000 3.7842 0.1697 .990 .980 0.6406 — 0.0943 
16-1200 6.6178 0.3483 986 .972 1253 1938 
17-0000 5.6933 0.1384 997 .994 .9578 0770 
17-1200 6.0083 0.2198 .993 .986 0154 1224 
18-0000 6.0189 0.1978 995 .990 0630 1100 
18-1200 5.1013 0.1940 993 .986 8650 1083 
19-0000 6.0881 0.1843 .995  .990 0284 1028 
19-1200 2.6928 0.2074 972  :945 4666 1158. 
20-0000 3.3309 0.2818 963 .927 5644 1593 
20-1200 2.2316 0.1834 965 .931 3754 1031 
21-0000 5.3035 0.2107 991 .982 8700 1187 
Ze 200 3.1120 0.3134 .953 .908 5487 1746 
22-0000 4.4330 0.1457 995 .990 .7500 0813 
22-1200 3.3728 0.1415 991 .982 5573 0802 
23-0000% 0.6804 0.3303 .614 .377 1954 1649 
23-1200% 0.8355 0.2358 746 .557 1885 1316 
TIME MEAN 4.0815 906 


STD. DEV. 1.8749 
mes A) 


*Corresponds to Ryo < .85. 





TABLE XILEL 


Volwessonwen ond related statistics 
for the data-sample for Spokane, Wash. 











STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE AX -VALUE 2 
DAY-GMT CALCULATED A -VALUE Rig = Ry a SY IX 
16-0000 6.1242 0.2171 994 .988 0350 1209 
16-1200 5.3894 0.0292 .999 .998 .9060 0163 
17-0000 6.2015 0.2579 992 .984 .0971 1432 
17-1200 5.4573 0.0661 .999 .998 .9586 0367 
18-0000 4.1476 0.2953 .973  .947 6969 1673 
18-1200 3.1086 0.2754 959  .920 5305 1562 
19-0000 3.3134 0.1454 .990 .980 5484 0825 
19-1200 2.5267 0.1350 .985  .970 4199 0766 
20-0000 2.6138 0.1618 .980 .960 4352 0914 
20-1200 4.5322 0.4601 948  .899 7780 2593 
21-0000 4.2112 0.3912 956 .914 .7166 2204 
P1-1200 3.6604 0.2030 .984 .968 6057 1145 
22-0000* 1.6059 0.3356 .822 .676 3184 1894 
22-1200% 1.0266 0.2085 .829 .687 2015 1176 
23-0000 5.3603 0.2307 .990 .980 8799 1299 
P1200 4.0725 0.4338 968 .937 5157 1406 © 
TIME MEAN 3.9594 925 


SED, DEV. 1.5414 
A 
“Corresponds to Ryo < .85. 





TABLE XIV 


Values of A and related statistics 
for the data-sample for Quillayute, Wash. 








STANDARD CORRELATION STANDARD STANDARD 
ERROR OF COEFFICIENT DEVIATION ERROR 
DATE NX -VALUE 2 
DAY-GMT CALCULATED A -VALUE Ria Ry ay Sy 1x 
16-0000 5.5314 1401 996 .992 0899 0.0992 
16-1200 5.5505 1507 996.992 1368 0.1065 
17-0000 4.3589 2506 979.958 8494 0.1789 
i200" 1.4652 3352 771.594 3626 0.2394 
18-0000% 1.7747 3777 .793-.629 4275 0.2701 
18-1200 1.6247 2586 867.752 3580 0.1850 
19-0000 2.0580 1481 .968  .937 4061 0.1058 
19-1200 2.3200 1796 963.927 4596 0.1282 
20-0000 2.5319 1793 969  .939 4977 pms 
20-1200 3.6523 3188 954  .910 een 0.2272 
21-0000* 1.6818 2995 842.709 3805 0.2133 
pe 200 1.4532 . 3086 joe 6a 3485 0.2199 
22-0000 1.7184 2594 Bi eye 0.1846 
22-1200 2.0152 1758 DE) O60 4858 0.1019 
23-0000 4.1114 1946 986 .972 7917 0.1382 
23-1200 3.9581 3402 972.945 7002 0.1755 
TIME MEAN 2.8629 851 


op, DEV. 1.5414 
7.C A ) 
*Corresponds to Rj2 < .85. 





faetable Til, a Yow correlation of 0.642 occurs at 22 March, OOQOOGMT, 
which results in a relatively high standard error Sy|x. Note in Table III 
that most of the profiles have Ry» >0.95 with associated small standard 
errors. There were two cases in Table III with low correlation coefficients, 
Specifically 0.704 and 0.642. This pattern of correlation coefficients 
R19 2 0.85 except for a possible set of one to three profiles for which 
R42 fell below this value, continues throughout all of the 12 stations 
considered here as shown in Tables III through XIV. Specific instances in 
these tables were the R,9-value value drops below 0.85 have been marked 
by an asterisk alongside of the date. In all of these cases the parameter 

AX is likewise smaller than normal in view of Eq. (20) 
h = Rio Fy / yx) (44) 

All of these special cases where R12 £0.85 are characterized by a 
Beamderd deviation % y lower than typical for the station. Such cases of 
swat *y-values represented cases where statistically there was insufficient 
vertical variation in the moisture range to ascribe it successfully to the 
power-profile parameterization. 

An example of a smaller-than-normal correlation coefficient Rio and 
its correspondingly small A has been extracted from the regression print 
out in the form of Table XV. This table lists the statistical results 
from Waycross, Ga., at OOOOGMT, 22 March, 1971. This power-profile 
analysis resulted in Ryo = 0.642 and <A = 0.9815, which was considerably 
lower than the mean A -value. In Table XV, the observed moisture profile 
is indicated in column 2 as Y of Eq. (7), and in column 3 is shown the 
Peeresponding ¥ On Hae 40): The residual shown in column 4 results from 
forming y-2 from Eq. (41). As expected, g began as a maximum positive 


value at the surface and decreased to 0.000 at 500 mb and becomes negative 
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TABLE XV 


An example for Waycross, Ga. (0000 GMT, 22 Mar., 1971) 
showing a condition leading tuo an ill-defined 4 -profile. 





PRESSURE (mb) Y (obs) 7 RESIDUAL 
1014 0.4660 0.3014 0.1646 
1000 0.3179 0.2955 0.0224 

950 Oe27 19 0.2736 =-0.0017 
900 0.2478 0.2505 -0.0028 
850 0.2020 0.2262 -0.0241 
800 0.1024 0.2003 -0.0979 
750 -0.0060 0.1728 -0.1789 
700 -0.1507 0.1434 -0.2942 
650 -0.1579 0.1118 =On27 le 
600 . 0.1296 Ona 0.0519 
550 - -0.0082 0.0406 -0.0488 
500 0.0000 0.0000 0.0000 
450 -0.3650 -0.0449 -0.3201 
400 -0.6949 -0.0951 -0.5998 
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TABLE XVI 


Environmental parameters used in the 
calculation of Table XV. 





RELATIVE MIXING 
PRESSURE (mb) TEMP (C) HUMIDITY, % RATIO (g/kg) 

1014 16.8 33 3.8801 
1000 1S eiG 25 2.7590 
950 220 2G 2.4816 
900 8.3 31 PSE Ty) 
850 6.9 29 2.1130 
800 iwa3 26 1.6800 
750 1.6 23 1.3087 
700 ieee 5 18 .9378 
650 - 2.9 20 .9186 
600 . - 6.3 48 1.7885 
550 -10.9 48 a soil 
500 -16.9 77 1.3270 
450 -19.7 39 .5726 
400 -27.2 34 2678 


rae Pp Pi > pa = TT 
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above 500 mb. The observed profile, Y(p), began also as a positive value 
at the surface, but becomes negative in the layer 750 to 600 mbs indicating 


that the actual mixing ratio profile had values less than w599 in that layer. 


Table XVI lists the actual atmospheric-sounding parameters which resulted 
in the values shown in Table XV. Note that the reported relative humidities 
are low in the 750 to 650 mb layer and high in the layer 600 to 450 mb. 
Since by Eq. (9), w is a function of relative humidity as well as of 
temperature and pressure, relatively low w-valuesexist from 750 to 650 mb 
and high w-valuesin the 600- to 450 mb- range. This, at least locally, is 
inconsistent with the assumption on which the lower-fit is based and, as 
Belated, gives a low value of Rj5. Note that the power-fit need not 
have been inconsistent if the temperature lapse was large enough through 
the 600 - 450 mb layer, which in this instance was not the case. 

As just noted, the second major factor which controls the consistency 
of the profile-fit at level p is the temperature T(p). It can also be 
Shown, by example, that a temperature inversion or an isothermal layer 
may be (depending upon relative humidity) associated with w values which 
do not decrease in height, resulting in the low correlation coefficients 
Rj 9- In Tables III to XIV, all cases where Rj, 2 0.85 which thus led to 
an ill-defined (i.e., too small) value of A are denoted by an asterisk. 
The explanation of all of these special cases is similar to that advanced 
for the Waycross, Ga. case in Tables XV, XVI. 

Note finally that it is convenient to list the time-mean A and Ry» 
fomeeedcin Station listed in Tables III,..., XIV. In addition it is desir- 
able to compute the temporal standard deviation of the A's in each 
E4b1e according to 


os 
; 
& 


nae) = im A, - A ech 16 = 
15 





(45) 
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B. REDUCTION OF VARIANCE OF A 


Having obtained the sixteen A -values per station the next portion 
of the investigation involved determining the feasibility of reducing the 
variance of <A for each station over the period of time considered using 
other sounding parameters. The synoptic variables introduced into the 
computation as independent variables permitted utilization of the stepwise 
multiple regression aspect of the BIMED 02R program where A _ was the 
dependent variable. In addition to the 16 values of XA per station, 16 
values of each of the following x. were read into the program as indepen- 
dent variables: 

surface temperature defined as xy 
surface mixing ratio defined as X» 
850 mb temperature defined as X, 
850 mb mixing ratio defined as X 


. 4 
700 mb temperature defined as X 


5 

500 mb temperature defined as X;. 
The above variables were used in the multiple regression as inputs with 
the computed A as the dependent variable in each case. The desired 
solution equation would then be of the form of Eq. (32) where Yy = A, 
and the Xx. are the respective independent variables. The stepwise 
regression program BIMED OZR was then utilized to derive the best-fit 
form of 


View), aaa a (46) 


0 ae 2 
Among the factors which might be considered as useful in reducing the 


temporal variance of XA are variations in the synoptic patterns, and those 


which may reflect purely local conditions. Among the latter effects 
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[meee tables JJ¥,..., XIV), station location and elevation should be 
considered, as well as proximity to large water bodies and lakes. In 

the sounding variables eae the lower-level parameters, were con- 
Sidered to be primarily representative of local effects, whereas the 
higher-level variables were presumed to be indicative of the larger-scale 
synoptic time-scales, although there is undoubtedly some interaction 
between the low- and high-level variables. 

The temporal standard deviations coat X ) listed at the base of 
Tables III,...,XIV resulting from Eq. (45) ranged from minimal values 
of oe A ) = 0.694 at Sault Ste Marie, and oa X ) = 0.701 at 
San Diego to a maximum 0,( A ) = 1.875 at Glasgow, Montana. In general 
it can be seen that Stations close to large bodies of water (Charleston, 
San Diego, Maniwaki, and Sault Ste Marie) have relatively small standard 
deviations whereas stations away from bodies of water have larger standard 
deviations (Bismarck, Glasgow, and Spokane). The cases of large values 
of a al X ) seem to afford considerable opportunity for utilization 
of the regression procedure for the reduction of the AX -variance. How- 
ever, it was decided to test the diagnostic usefulness of Eq. (46) for all 
twelve stations, regardless of the aed X +) value. 

A preliminary test run was conducted using the BIMED 02R program with 
the variables defined previously in the section. During this run all 
independent variables were allowed to enter freely based solely on their 
partial correlation values. The program was allowed to run until the 
requirement of 2% added fractional explained variance after introduction 
of the kth variable could not be met. 

Two facts became evident from this preliminary test: (1) the indepen- 


dent variable with the highest initial correlation coefficient in almost 
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all cases was X,, and (2) the maximum number of significant steps necessary 
to ensure the required minimum added explained variance at each step was 
four. At steps five and six of the preliminary run the added explained 
variance upon predictor-entry was too small. Since a low-level moisture 
predictor was desired to be compatible with the variations of the A - 
profile parameter (cf., the discussion for Waycross at 0000 GMT, 22 March, 
1971), it was decided to utilize the control-delete option of BIMED 02R 
Semerce in A, aS Che first predictor at each station. The second option 
used was that of limiting the number of predictors to four in each 
regression analysis. The results of the sets of regressions over the 12 
stations are shown in Table XVII. 

Equation (42) gives the relationship between the standard error of 

nN , after the four-step multiple correlation R. has been determined 
S ayy GRE OCA) (47) 

From Table XVII it can be seen that when oo XA ) is relatively small, 
the explained fraction of this statistic is relatively small. The result, 
then, of application of the multiple regression technique is to reduce 
the standard deviation kk X +) of all stations, but to reduce the 
Stations with the largest Are X ) the most. 

San Diego is a good example of a station having a small or eee 

es ( XA ) = 0.7012, which means, according to the above result, a 

relatively small reduction, specifically to S ypx=0.5686. Correspondingly, 
this small additional reduction is reflected in the relatively low multiple 
Correlation coefficient at San Diego of eS .5852. At the other extreme 
is Glasgow having a temporal standard deviation ise XA ) = 1.8749. By 
the multiple regression procedure a standard error of estimate of 0.8745 
results. Again this is consistent with the high multiple correlation 
Poeretcrent, KR, = 0.8846. 
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OF STATIONS 


NUMBER 


.30 3 5 .60 -65 70 75 .80 85 .90 95 


STANDARD ERROR 


Figure 2: Histogram depicting number of stations within certain 
standard error limits. 


ai, 





miececesult or: the multiple regression for each station is to give a re- 
duced standard error of A. A histogram depicting the 12 standard error 
Values is shown in Figure 2, on a scale of S AL X Prom 0750 to 0.95. Note 
Giememost of the Standard errors after regression fall into the range 
eto 0.65. 

The righthand column in Table XVII lists the order of entry of the 
independent variables into the stepwise regression and is based on the 
value of the partial correlation of the respective variables as explained 
in Section III (B). From Table XVII it may be noted that in seven of the 
12 stations the second independent variable added to the regression was Xo. 
Also in seven of the 12 stations the third independent variable added was 
Poeeetmas led to the idea of possible classification of stations according 
to a forced order of entry, an option of BIMED O02R, and will be discussed 


in the next sub-Section. 


C. STATION CLASSIFICATION ACCORDING TO THE PREDICTOR ENTRY 


ihe AVIIT lists an attempt to classify stations by the relative 
priority of the independent variables to the program specified in the order 
Xn, Xe, Xy > and the remaining predictor in a 4-step selection left as free 
choice according to the selection criterion. As can be seen by comparing 
Tables XVII and XVIII there was no change in the multiple correlation 
coefficient at 10 of the 12 stations as a result of specifying the 4-6-2 
fiepure order. This 15 because in almost alll cases the variables X,, X¢, and 
X9 were in the original run (Table XVII) even though not necessarily in this 
specific order. Table XVIII lists the "single free choice'' variable under 


order of input as the last independent variable in the list for each station. 


Dll 
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It may be noted that there are four stations with surface temperature as 
the last independent variable (Charleston, Maniwaki, Bismarck, and 
Spokane), four with 500 mb-temperature (Waycross, Lake Charles El Paso, 
and Caribou), and four with 850 mb-temperature (San Diego, Sault Ste 
Marie, Glasgow, And Quillayute). There might, then, be some way of 
classifying or typing stations on the basis of the free-choice predictor 
selected at step 4. However, with the small station set available it 
was not clear that this could be done. 

Finally, it is to be emphasized that Table XVII ree that the 
observed A -values tend to adjust to (a) the synoptic-scale time 
variations in the sounding and to (b) the local-effect variations, The 
Mieeer eitect is harder to specify so that one achieves some limiting 
accuracy in the expected value resulting from the regression. The magnitude 
of this accuracy is suggested by the standard error which tends to show a 


cluster of values in the range 0.50 to 0.65 (Fig. 2). 


SN) 





Venere tiAbLE waleR VAPOR TESTS 


As noted in the introduction, two verification techniques were to be 
M@sed Co evaluate the results of this investigation. The first of these, 
namely verification of the moisture profile, has been discussed in detail 
in Section IV from a statistical viewpoint. In this section, an analysis 
of the verification of the moisture-parameterization technique from the 
viewpoint of precipitable water vapor is presented. 

The formulas used to calculate the precipitable water are shown as 
Eqs. (12) and (13). Equation 12 was utilized to obtain precipitable water 
based upon the moisture profile of the data for the given data sample. It 
is the precipitable water from the observed sounding, and hence will be 


referred to as observed precipitable water, W Equation 13 is used to 


Ae 
Sebeulate Che precipitable water utilizing the A -values already obtained. 
As can be seen from Eq. (13), it is based on only three variables: 500-mb 
mixing-ratio, surface pressure, and A. There are, however, two possible 
methods since A -values were obtained by the profile~parameter method 
and, secondly, by the Poetic PetmesctOne tecimague.  Préecipicable water 
calculated using the profile parameter <A -values, one for each sounding, 
will be Ree oabed as W., and those relatively smooth values obtained by 
using the multiple regression 4 -values denoted by W... 

For this phase of the investigation, four stations were chosen primarily 
on the basis of their varying geographical locations: Waycross, Ga., in 


the southeastern United States; San Diego in the southwest; Spokane, Wash., 


in the northwest; and Sault Ste Marie, Mich., in the northeast. Results 
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@iethne precipitable water calculations are listed in Tables XIX through 
Peteand are plotted in Figures 3 through 6. 

As can be seen from Tables XIX,...,XXII or from Figures 3,...,6, the 
Meeevalue is in almost all cases within approximately 0.05 cms of the W, 
value. Needless to say, with such good verification, the trend of Wee 
is the same as the observed trend. It is obvious from the results, then, 
Phat the use of the power-profile technique can be expected to yield 
precipitable water values which are very representative of the W, values, 
meeevrthin 0.05 cms (cf., Tables XIX,...,XXII). 

The values of Wop are also presented in Tables XIX,...,XXII and in 
Figures 3,...,6. As can be seen errors do exist in the W,, values, 
assuming W, is the correct value. These errors can at times become large, 
i.e. as much as 100%. In general, however, the W., values which result 
from using A -values obtained by multiple regression follow the time 
trend of Noe Large errors in We, are then the exception rather than the 
rule, and as can be seen in Figures 3,...,6, occur at only one or two 
sounding times per station. The soundings which give the largest Woy 
errors are, in fact, just those soundings which had the lower values of 

r and of the correlation coefficients R,j2 in connection with the Wes, 
computations. 

The 0000 GMT, 22 March, 1971 sounding at Waycross, Ga. will be used 
to illustrate the point of a large Wem error. As can be seen in Table IIIT, 
this sounding gave a low correlation coefficient value, 0.642, anda A - 
watue ot 0.9815 where the timeé-mean over the 16 soundings for that station 
was 2.8557. Figure 3 shows that the use of the A -value from the power- 
fit profile, A = 0.9815, in Eq. (13) gives an accurate verification of 


Wo. On the other hand the value of A resulting from the multiple 
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TABLE XIX 





Table of precipitable water vapor values (Wo Maas. the), 
for Waycross, Ga. 
DAY OBSERVED PROFILE MULT. REG. 
DAY - GMT Ne ee Mee 





16 - 0000 7190 3.8704 0169 
16 - 1200 0227 1884 2003 
17 - 0000 9360 8618 4918 
17 - 1200 5325 Sal 5940 
18 - 0000 5678 5863 5923 
18 - 1200 6549 6600 5514 
19 - 0000 1639 Rilava RA 
19 - 1200 2758 2377 9609 
20 - 0000 One 9818 77 
20 - 1200 5083 5216 5058 
21 - 0000 4897 PRE 5948 
pi 1200 8447 8370 9129 
22 - 0000 9807 1674 4002 
22 - 1200 2914 A747 fess 
23 - 0000 5845 3996 fo 
23 - 1200 0389 8514 1848 





TABLE XX 


Table of precipitable water vapor values (Wo, Wes, Wem) 
LOEeoamewneoO. "Ga lit. 





DAY OBSERVED PROFILE MULT. REG. 
DAY - GMI Wo Wes Wom 
16 - 0000 . 7488 . 7849 5459 
16 - 1200 .8116 .8317 6233 
17 - 0000 ,8222 . 8833 6292 
17 - 1200 . 8586 . 8796 .8950 
18 - 0000 .8783 8523 .7305 
18 - 1200 OSE 0334 .7619 
19 - 0000 .9586 rosy . 8904 
19 - 1200 .1637 .2238 6504 
20 - 0000 .7991 .9855 .5071 
20 - 1200 .7310 . 8049 .8771 
21 - 0000 - 6700 £7594 .8141 
21 - 1200 9416 .9233 . 8065 
22 - 0000 205 8268 .0792 
22 - 1200 .2398 3404 . 8782 
23 - 0000 .1738 3048 .7369 
23 - 1200 7495 . 7309 4738 


Ne A a ae 
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TABLE XXI 


Table of precipitable water vapor values (Wo, Wes, Wen) 
for Spokane, Wash. 


DAY 
DAY - GMT 
16 0000 
16 1200 
1g 0000 
17 1200 
18 - 0000 
18 1200 
lig 0000 
ng 1200 
20 0000 
20 1200 
Zh 0000 
21 - 1200 
22 - 0000 
jag 1200 
Ze 0000 
5 1200 


en PP SS 


OBSERVED 


Wo 


10970 
. 3309 
70. 
. 3807 


3491 


~4091 


-4692 


eo" 


eae) 


.6018 


-4113 


. 3094 


so076 


O72 


. 1988 


22039 


58 


PROFILE 
Wes 


By be ho 4 
(33753 
~4002 
MEOH 
~3174 
1377-0 
4820 
9436 
ye So 
-4964 
3489 
2940 
- 3036 
wee 7 
foo 5 


.0635 


MULT. REG. 
Wom 


mee 2 


a2 70 


~4135 


nOO97 


egos 


OD 


4334 


OG) 


7OOS2 


3805 
3469 
. L8a8 
~4224 
5 Shera: 
6473 


13730 





TABLE XXI 


Table of precipitable water vapor values (Wo, Wes, Wen) 
for Sault Ste Marie, Mich. 





DAY OBSERVED PROFILE MULT. REG. 
DAY - GMT Wo Wes Wem 
16 - 0000 . 1656 1.1476 .0986 
16 - 1200 sBi12 3207 3854 
17 - 0000 13920 3293 PS 
17 - 1200 3501 .3401 .3188 
18 - 0000 3865 3693 .3514 
18 - 1200 .4089 . 3887 2807 
19 - 0000 .5890 .5605 . 7646 
19 - 1200 . 7004 .7405 . 8876 
20 - 0000 .7978 .7349 6034 
20 - 1200 po nlole .5162 4981 
21 - 0000 5046 SES, .4953 
2 = 1200 2910 2889 .4080 
22 - 0000 Bie 3433 .4346 
22 = 1200 252 O00 2942 
23 - 0000 3232 75234 poz 
23 - 1200 .2074 #2032 .1921 
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regression procedure for the date and station considered was A = 2.6793. 
Since the latter value is considerably larger than the sounding best-fit 
value of A, the application of the multiple regression , -value in 

Eq. (13) gives a much larger value of Wom than Wo? hence the large error, 


Slee - Wo; 


at this sounding time (Figure 3). 

The results of Tables XIX,...,XXII (and the four supporting diagrams 
Figures 3,...,6) have been presented in the statistically summarized form 
as Table XXIII. In Table XXIII, the mean W,7values have been compared 
to the computed mean values of W., and W. using the 16 sounding samples 
at the four stations, Waycross, San Diego, Spokane, Sault Ste Marie. 


The results were by-products of the correlation of W, separately against 


W and of Wo with W 


cs? which was applied to the two sets of data in 


cm? 
Maples XIX,...,XX11, respectively. 

The standard deviations of Wo, Wes, and Wem have been summarized in 
Table XXIII (column 3). Note that the means and standard deviations of 
the W-samples computed by the different procedures are nearly equal. 

ie correlation coefficients R(W,, Was) and R(W,, W,,,) are listed for 
the four Stations in column 4 of Table XXIII. The corresponding 
coefficients of determination R* are listed in column a, and aves cto. be 
interpreted as the fractional explained variance (of Wo): Finally the 


Pomumi.otd. Error’ is to be interpreted as the mean dispersion of the 


curve of Wo minus W 


cg» and also that of W, minus W.., using the following 


4 


equation for standard error: 


2 
SWo|We = CL - Re FW, (48) 
Note that the estimation of Weg to give W, leads to a standard error (a 


mean residual) of 0.07 cm, whereas the corresponding standard error in the 
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comparison of Wes with W, is close to 0.31 cms. This larger error arises 
Medmarily from the southerly stations in Table XXIII. 

In conclusion, the value of Wes based upon a best-fit AX -value in 
conjunction with Eq. (13) gives a good estimation of ee On the other hand, 
the computation of W.,, wherein the A -value is estimated by a multiple 
regression formula uSing gross parameters from the sounding has some skill 
in estimation of ee but perhaps needs additional study for selection of 
input variables for optimal inclusion in the multiple regression equation 


used to determine A. 
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VI. CONCLUSIONS 


At a group of 12 geographically diverse stations, the mixing-ratio 
profile was parameterized as a power-law profile, Eq. (6), having the 
characteristic exponent A . For each sounding A _ was derived by a 
least-squares fit with an accuracy usually characterized by a correlation 
coefficient of 0.95 or higher. At each station, the fit was most accurate 
in this form when the standard deviation of observed mixing-ratio was 
large. The set of A -values were grouped together in time series at each 
station, and they were found to be quite variable in general. This 
variability of A was related to the detailed structure of the moisture 
profile from the sounding data. A sizeable part of the standard devia- 
tion ore X ) was found to be explainable by use of a multiple regression 
procedure which incorporated easily defined properties of the sounding. 
The primary test of the A -profile parameterization was made in application 
to the comparison of observed versus computed values of the precipitable 
water vapor at four stations, with the latter given by 


(W500) (500) | /Ps \ A+1l [Ptop \ Atl 
Seen + 1) seo 500 (13) 


On an individual sounding basis the X -value gave good verification 
of the observed precipitable water vapor value. Over a period of time 
it was possible to generate multiple regression equations for XA which 
gave moderately good verification when applied for the appropriate sound- 
ing using Eq. (13). There is some possibility that this procedure can be 
employed in a predictive sense provided the parameters used in the 


regression can be successfully generated by a primitive equations model. 
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It is possible that the multiple regression sounding parameters may be 
selected to yield closer agreement with the j -values determined from 


the individual soundings. 
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